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A solution to the more than 300-years old problem of geometric and 
physical interpretation of fractional integration and differentiation (i.e., 
integration and differentiation of an arbitrary real order) is suggested 
for the Riemann-Liouville fractional integration and differentiation, the 
Caputo fractional differentiation, the Riesz potential, and the Feller 
potential. It is also generalized for giving a new geometric and physical 
interpretation of more general convolution integrals of the Volterra 
| type. Besides this, a new physical interpretation is suggested for the 

Stieltjes integral. 
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1 Introduction 

X 

It is generally known that integer-order derivatives and integrals have clear 
physical and geometric interpretations, which significantly simplify their use 
for solving applied problems in various fields of science. 

However, in case of fractional-order integration and differentiation, which 
represent a rapidly growing field both in theory and in applications to real- 
world problems, it is not so. Since the appearance of the idea of differentia- 
tion and integration of arbitrary (not necessary integer) order there was not 
any acceptable geometric and physical interpretation of these operations 
for more than 300 years. The lack of these interpretations has been ac- 
knowledged at the first international conference on the fractional calculus in 
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New Haven (USA) in 1974 by including it in the list of open problems [21]. 
The question was unanswered, and therefore repeated at the subsequent 
conferences at the University of Strathclyde (UK) in 1984 [15] and at the 
Nihon University (Tokyo, Japan) in 1989 [19]. The round-table discussion 
[13, 10, 14] at the conference on transform methods and special functions 
in Varna (1996) showed that the problem was still unsolved, and since that 
time the situation, in fact, still did not change. 

Fractional integration and fractional differentiation are generalisations 
of notions of integer-order integration and differentiation, and include n-th 
derivatives and n-folded integrals (n denotes an integer number) as partic- 
ular cases. Because of this, it would be ideal to have such physical and 
geometric interpretations of fractional-order operators, which will provide 
also a link to known classical interpretations of integer-order differentiation 
and integration. 

Since the need for the aforementioned geometric and physical interpre- 
tations is generally recognised, several authors attempted to provide them. 
Probably due mostly to linguistical reasons, much effort have been devoted 
to trying to relate fractional integrals and derivatives, on one side, and frac- 
tal geometry, on the other [18, 27, 9, 16, and others]. However, it has been 
clearly shown by R. Rutman [22, 23] that this approach is inconsistent. 

Besides those "fractal-oriented" attempts, some considerations regarding 
interpretation of fractional integration and fractional differentiation were 
presented in [16]. However, those considerations are, in fact, only a small 
collection of selected examples of applications of fractional calculus, in which 
hereditary effects and self-similarity are typical for the objects modelled with 
the help of fractional calculus. Although each particular problem, to which 
fractional derivatives or/and fractional integrals have been applied, can be 
considered as a certain illustration of their meaning, the paper [16] cannot 
be considered as a definite answer to the posed question. 

A different approach to geometric interpretation of fractional integration 
and fractional differentiation, based on the idea of the contact of a-th order, 
has been suggested by F. Ben Adda [1, 2]. However, it is difficult to speak 
about an acceptable geometric interpretation if one cannot see any picture 
there. 

Obviously, there is still a lack of geometric and physical interpretation 
of fractional integration and differentiation, which is comparable with the 
simple interpretations of their integer-order counterparts. 

In this paper we present a new approach to solution of this challenging 
old problem. 

We start with introducing a simple and really geometric interpretation 
of several types of fractional-order integration: the left-sided and the right- 
sided Riemann-Liouville fractional integration, the Riesz potential, and the 
Feller potential. 

Based on this, a physical interpretation of the Riemann-Liouville frac- 
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tional integration is proposed in terms of inhomogeneous and changing (non- 
static, dynamic) time scale. Moreover, on this way we give a new physical 
interpretation of the Stieltjes integral. We also try to persuade the readers 
that the suggested physical interpretation of fractional integration is in line 
with the current views on space-time in physics. We also suggest physical 
interpretation for the Riemann-Liouville fractional differentiation and for 
the Caputo fractional differentiation. Finally, we show that the suggested 
approach to geometric interpretation of fractional integration can be used 
for providing a new geometric and physical interpretation for convolution 
integrals of the Volterra type. 

2 Geometric interpretation of fractional 
integration: Shadows on the walls 

In this section we first give a geometric interpretation of left-sided and right- 
sided Riemann-Liouville fractional integrals, and then consider the Riesz 
potential. 

2.1 Left-sided Riemann-Liouville fractional integral 

Let us consider the left-sided Riemann-Liouville fractional integral [20, 24] 
of order a, 

t 

o^f(t) = f ^- ) Jf(r)(t-rr- 1 dr, (1) 
and write it in the form 

f 

/f/(t) = J f(r)dg t (r), (2) 
o 

The function gt(r) has an interesting scaling property. Indeed, if we take 
t\ = kt and t\ = kr, then 

9t 1 (n)=9kt(kT) = k a g t (r). (4) 

Now let us consider the integral (2) for a fixed t. Then it becomes simply 
a Stieltjes integral, and we can utilize G. L. Bullock's idea [3]. 

Let us take the axes r, g, and /. In the plane (r, g) we plot the function 
gt{j) for < r < t. Along the obtained curve we "build a fence" of the 
varying height /(r), so the top edge of the "fence" is a three-dimensional 
line (T,g t (T),f(T)),0<T<t. 
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This "fence" can be projected onto two surfaces (see Fig. 1): 

• the area of the projection of this "fence" onto the plane (r, /) corre- 
sponds to the value of the integral 

t 

o4(t) = / /(r)dr; (5) 
o 

• the area of the projection of the same "fence" onto the plane (g, f) 
corresponds to the value of the integral (2), or, what is the same, to 
the value of the fractional integral (1). 

In other words, our "fence" throws two shadows on two walls. The first 
of them, that on the wall (r, /), is the well-known "area under the curve 
/(t)", which is a standard geometric interpretation of the integral (5). The 
"shadow" on the wall (g, f) is a geometric interpretation of the fractional 
integral (1) for a fixed t. 

Obviously, for gt(r) = r both "shadows" are equal. This shows that 
classical definite integration is a particular case of the left-sided Riemann- 
Liouville fractional integration even from the geometric point of view. 

What happens when t is changing (namely growing)? As t changes, the 
"fence" changes simultaneously. Its length and, in a certain sense, its shape 
changes. For illustration, see Fig. 2. If we follow the change of the "shadow" 
on the wall (g,f), which is changing simultaneously with the "fence" (see 
Fig. 3), then we have a dynamical geometric interpretation of the fractional 
integral (1) as a function of t. 

2.2 Right-sided Riemann— Liouville fractional integral 

Let us consider the right-sided Riemann-Liouville fractional integral [20, 24], 

6 

and write it in the form 

b 

t lZf(t) = J f(r)dh t (r), (7) 
t 

Then we can provide a geometric interpretation similar to the geometric 
interpretation of the left-sided Riemann-Liouville fractional integral. How- 
ever, in this case there is no any fixed point in the "fence" base - the end, 
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Figure 2: The process of change of the fence basis shape 
for It /(*), a = 0.75. 
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Figure 3: Snapshots of the changing "shadow" of changing 
"fence" for I?f(t), a = 0.75, f{t) = t + 0.5sin(t), 
with the time interval At = 0.5 between the snap- 
shops. 



corresponding to r = 6, moves along the line r = b in the plane (r, g) when 
the "fence" changes its shape. This movement can be observed in Fig. 4. 
(In the case of the left-sided integral, the left end, corresponding to r = 0, 
is fixed and does not move.) 

All other parts of the geometric interpretation remain the same: the 
"fence" changes its shape as t changes from to 6, and the changing shad- 
ows of this "fence" on the walls (g, f) and (r, /) represent correspondingly 
the right-sided Riemann-Liouville fractional integral (6) and the classical 
integral with the moving lower limit: 

6 

t4(t) = J f(r)dr; (9) 
t 

Obviously, for <7t(r) = r both "shadows" are equal. Therefore, we see 
that not only the left-sided, but also the right-sided Riemann-Liouville frac- 
tional integration includes the classical definite integration as a particular 
case even from the geometrical point of view. 
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Figure 4: The process of change of the fence basis shape 
for t /f /(t), a = 0.75. 



2.3 Riesz potential 

The Riesz potential [20, 24] 





oRZf(t) = f ^- ) Jf(r)\T-tr 1 dT 



(10) 



is the sum of the left-sided and the right-sided Riemann-Liouville fractional 
integrals: 



t 



o 



r(a) 



f(r)(r-tr- l dr. (11) 



o t 
The Riesz potential (10) can be written in the form 

b 



oRbf(t) = / f(r)dr t (r), 



o 



n(T) = T(a\l) { ta + Sign(T " t] |T " * r }' 



(12) 



(13) 



The shape of the "fence", corresponding to the Riesz potential, is de- 
scribed by the function tj(t). In this case the "fence" consists of the two 
parts: one of them (for < r < t) is the same as in the case of the left-sided 
Riemann-Liouville fractional integral, and the second (for t < r < 6) is the 
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Figure 5: The process of change of the fence basis shape 
for the Riesz potential oRf f(t), a = 0.75. 



same as for the right-sided Riemann-Liouville integral, as shown in Fig. 5. 
Both parts are joined smothly at the inflection point r = t. 

The shape of the "fence" , corresponding to the Riesz potential, is shown 
in some of its intermediate position by the bold line in Fig. 5. Obviously, 
Fig. 5 can be obtained by laying Fig. 4 over Fig. 2, which is a geometric 
interpretation of the relationship (11). 

The shadow of this "fence" on the wall (g, f) represents the Riesz po- 
tential (10), while the shadow on the wall (r, /) corresponds to the classical 
integral 



For a = 1 both "shadows" are equal. This shows that the classical definite 
integral (14) is a particular case of the Riesz fractional potential (10) even 
from the geometric point of view. We have already seen this inclusion in the 
case of the left-sided and the right-sided Riemann-Liouville fractional inte- 
gration. This demonstrates the strength of the suggested geometric inter- 
pretation of these three types of generalization of the notion of integration. 

2.4 Feller potential 

The Feller potential operator & a f(t) is, similarly to the Riesz potential, also 
a linear combination of the left- and right-sided Riemann-Liouville fractional 
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(14) 



o 
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integrals, but with general constant coefficient; 



is c,d [24, Chap. 3]: 



$ a f(t) = c a I?f(t)+d t lZf(t). 



(15) 



The geometric interpretation of the Feller potential can be easily ob- 
tained by properly scaling and then superimposing Fig. 4 and Fig. 2. The 
"fence" obtained in this way is, in general, discontinuous at r = t. Its 
shadow on the wall (r, /) is equal to the classical definite integral (14). The 
shadow on the wall (g, f) consists, in general, of the two areas, which may 
overlap depending on the values of the coefficients c and d. 

3 Two kinds of time — I 

The geometric interpretation of fractional integration, given in the previous 
sections, is substantially based on adding the third dimension (for <7t(r)) 
to the classical pair r, f(r). If we consider r as time, then g(r) can be 
interpreted just as a "deformed" time scale. What could be the meaning of 
having - and using - two time axes? To answer this question, let us recall 
some facts of the history of the development of the notion of time. 

That were contributions of Barrows and Newton to the development of 
mathematics and physics in the XVII century which led to the appearance of 
the "mathematical time", which is postulated to "flow equably" and which 
is usually depicted as a semi- infinite straight line [26]. 

Newton himself postulated [17]: 

"Absolute, true and mathematical time of itself, and from its own 
nature, flows equably without relation to anything external." 

Such a postulate was absolutely necessary for developing Newton's dif- 
ferential calculus and applying it to problems of mechanics [26]: 

"The outstanding mathematical achievement associated with the 
geometrization of time was, of course, the invention of the cal- 
culus of fluxions by Newton." 

"Mathematically, Newton seems to have found support for 
his belief in absolute time by the need, in principle, for an ideal 
rate-measurer." 

The invention of differential and integral calculus and today's use of them 
is the strongest reason for continuing using homogeneous equably flowing 
time. 

Time is often depicted using the time axis, and the geometrically equal 
intervals of the time axis are considered as corresponding to equal time 
intervals (Fig. 6). 
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Figure 6: Homogeneous time axes. 
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Figure 7: Time slowing down. 



This assumption, however, cannot be neither proved nor rejected by 
experiment. Two lengths of geometric intervals can be measured and com- 
pared, since they are available for measurement simultaneously, at the same 
time (or, more precisely, at the same time and at the same place). Two 
time intervals can never be compared, because they are available to us for 
measurement (or for observation) only sequentially. 

Indeed, how do we measure time intervals? Only by observing some 
processes, which we consider as regularly repeated. G. Clemence wrote [7]: 

"The measurement of time is essentially a process of counting. 
Any recurring phenomenon whatever, the occurences of which 
can be counted, is in fact a measure of time." 

Clocks, including atomic clocks, repeat their "ticks", and we simply 
count those ticks, calling them hours, minutes, seconds, milliseconds, etc. 
But we are not able to verify if the absolute time which elapsed between, 
say, the fifth and the sixth tick (the sixth "second") is exactly the same as 
the time, which elapsed between the sixth and the seventh tick (the seventh 
"second"). This possible inhomogeneity of the time scale is illustrated in 
Fig. 7. 

The fact that time measurement as a process of counting of repeating 
discrete events does not really exclude inhomogeneity of time, has been 
nicely mentioned by L. Carroll in Alice's Adventures in Wonderland [6, 
Chap. 7]: 

"... I know I have to beat time when I learn music." 

"Ah! That accounts for it," said the Hatter. "He [Time] 
won't stand beating. Now, if you only kept on good terms with 
him, he'd do almost anything you liked to do with the clock. . . " 

Figures 6 and 7 show those "clock ticks", which we can register, only 
symbolically. One can interpret them as if there exists some absolute, or 
cosmic, inhomogeneous time axis, to which we can compare individual ho- 
mogeneous time, represented by some "clock ticks". Our picture of the 
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Table 1: Recording speed using slowing-down clocks 



Person N 


Recorded values 


Observer 


individual 


of velocity [m/s] 


absolute (cosmic) 


"seconds" 




"seconds" 





10 





1 


11 


1 


2 


12 


3 


3 


13 


7 


4 


12 


15 


5 


11 


31 


6 


10 


63 


7 


9 


127 



individual homogeneous time has the form shown in Fig. 6. The cosmic 
time may be not necessarily flowing equably, like that shown in Fig. 7. 

To illustrate the idea, let us consider the following situation. Suppose 
person N has two devices: one is a speedometer, and another one is the 
clock, which is slowing down, so the interval between the two subsequent 
ticks is double comparing to the interval between the previous ticks (see 
Fig. 7). Person N reads the velocity values indicated by the speedometer 
at each encountered "second", without knowing that the clock is, in fact, 
slowing down. 

Using these two series of data, namely the recorded sequence of values 
of speed, and the sequence of the counted "seconds", person N can estimate 
the distance which he has passed. 

For simplicity, let us suppose that the first "second" of the time shown by 
the clocks is equal to the absolute time "second" . The results of observations 
in this hypothetical experiment are given in Table 1. 

Person N will compute the distance he has passed as 

S N = 10 ■ 1 + 11 • 1 + 12 ■ 1 + 13 • 1 + 12 • 1 + 11 • 1 + 10 • 1 = 79. 

However, if there would be an independent observer O, knowing about 
the slowing-down clock, then such an observer would obtain a notably dif- 
ferent result for the distance passed by person N: 

S = 10 • 1 + 11 • 2 + 12 • 4 + 13 • 8 + 12 • 16 + 11 • 32 + 10 • 64 = 1368. 

Below we use this idea for giving a new mechanical interpretation of the 
Stieltjes integral. 
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4 Physical interpretation of 
the Stieltjes integral 

Imagine a car equipped with two devices for measurements: the speedometer 
recording the velocity v(r), and the clock which should show the time r. 
The clock, however, shows the time incorrectly; let us suppose that the 
relationship between the wrong time r, which is shown by the clock and 
which the driver considers as the correct time, on one hand, and the true 
time T, on the other, is described by the function T = g(r). This means 
that where the driver "measures" the time interval dr, the real time interval 
is given by dT = dg(r). 

The driver A, who do not know about wrong operation of the clock, will 
compute the passed distance as the classical integral: 

t 

S A (t) = J v(r)dr . (16) 
o 

However, the observer O knowing about the wrong clock and having the 
function g(r), which restores the correct values of time from the driver's 
wrong time t, will compute the really passed distance as 

So(t) = J v{r)dgir). (17) 
o 

This example shows that the Stieltjes integral (17) can be interpreted 
as the real distance passed by a moving object, for which we have recorded 
correct values of speed and incorrect values of time; the relationship between 
the wrongly recorded time r and the correct time T is given by a known 
function T = g(r). 

5 Physical interpretation of fractional integration: 
Shadows of the past 

Now let us consider the left-sided Riemann-Liouville fractional integral 

t 

So® = J v(r)dg t (r) = I?v(t), (18) 
o 

where gt(r) is given by (3). 

The fractional integral So(t) of the function v{t) can be interpreted as 
the real distance passed by a moving object, for which we have recorded the 
local values of its speed v (r) (individual speed) and the local values of its 
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time r (individual time); the relationship between the locally recorded time 
r (which is considered as flowing equably) and the cosmic time (which flows 
non-equably) is given by a known function ^(r). 

The function gt{j) describes the inhomogeneous time scale, which de- 
pends not only on r, but also on the parameter t representing the last mea- 
sured value of the individual time of the moving object. When t changes, the 
entire preceding cosmic time interval changes as well. This is in agreement 
with the current views in physics. Indeed, B. N. Ivanov [12, p. 33] mentioned 
that time intervals depend on gravitational fields. Similarly, S. Hawking [11, 
p. 32-33] wrote that: 

". . . time should appear to run slower near a massive body like 
the earth." [. . . ] 

". . . there is no unique absolute time, but instead each in- 
dividual has his own personal measure of time that depends on 
where he is and how he is moving." 

When a moving body changes its position in space—time, the gravita- 
tional field in the entire space-time also changes due to this movement. As 
a consequence, the cosmic time interval, which corresponds to the history 
of the movement of the moving object, changes. This affects the calculation 
(using formula (18)) of the real distance So(t) passed by such a moving 
object. 

In other words, the left-sided Riemann-Liouville fractional integral of the 
individual speed v(t) of a moving object, for which the relationship between 
its individual time r and the cosmic time T at each individual time instance 
t is given by the known function T = fff(r) described by the equation (3), 
represents the real distance So(t) passed by that object. 

6 Physical interpretation of the Riemann-Liouville 
fractional derivative 

On the other hand, we can use the properties of fractional differentiation and 
integration [20, 24] and express v(t) from the equation (18) as a left-sided 
Riemann-Liouville fractional derivative of So(t): 



where oDf denotes the Riemann-Liouville fractional derivative [20, 24], 
which is for < a < 1 defined by 



v(t) 



vD?S {t) 



(19) 




(20) 
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This shows that the left-sided Riemann-Liouville fractional derivative of 
the real distance So(t) passed by a moving object, for which the relationship 
between its individual time r and the cosmic time T at each individual time 
instance t is given by the known function T = gt(j) described by equation 
(3), is equal to the individual speed v{t) of that object. 

On the other hand, we can differentiate the relationship (18) with respect 
to the cosmic time variable t, which gives the relationship between the ve- 
locity vo(t) = S' (t) of the movement from the viewpoint of the independent 
observer O and the individual velocity v(t): 

v {t) = j t vl?v{t)= Dl- a v(t), (21) 

Therefore, the (1 — a)-th-order Riemann-Liouville derivative of the in- 
dividual velocity v(t) is equal to the velocity vo(t) from the viewpoint of the 
independent observer, if the individual time r and the cosmic time T are re- 
lated by the function T = gt(r) described by equation (3). For a = 1, when 
there is no dynamic deformation of the time scale, both velocities coincide: 
v (t)=v(t). 

7 Physical interpretation of the Caputo 
fractional derivative 

Applying fractional integration of order (3 = 1 — a to both parts of the 
relationship (21) gives: 

v(t) = O lt a vo(t) = l}- a S> o (t) = °D?So(t), (22) 

where ^Df denotes the Caputo fractional derivative [4, 5, 20], which is for 
< a < 1 defined by 

C oD?f(t) = -±—]P^. (23) 

tJy ' r(i-o) J (t-T) a v ' 

The relationship (22) is similar to (19). Therefore, the Caputo fractional 
derivative has the same physical interpretation as the Riemann-Liouville 
fractional derivative (see Section 6). This coincidence becomes more obvious, 
if we recall [20] that if /(0) = 0, then the Riemann-Liouville derivative 
and the Caputo derivative of order a (0 < a < 1), coincide: $Dff{t) = 
oD?f(t). 

8 Two kinds of time — II 

The suggested physical interpretation of fractional integration and fractional 
differentiation is based on using two kinds of time: the cosmic time and the 
individual time. 
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As mentioned above, due to the history of the development of mathemat- 
ics and physics, we are taught to think about the time, in fact, geometrically. 
The real roots of this go even far back to Euclid [26]: 

"Euclid considered space as the primary concept of science and 
relegated time to poor second." 

The entire integral and differential calculus is based on using mathemat- 
ical (homogeneous, equably flowing) time. There is no chance to change 
this state, and there is nothing to suggest instead of the classical calculus. 
Moreover, there is probably even no need for this. We can just realize that 
the classical calculus provides tools for describing the dynamic properties of 
the cosmic time, which - according to physicists - is inhomogeneous (flowing 
non-equably) . Indeed [11, p. 33-34], 

"The old idea of an essentually unchanging universe that could 
have existed, and could continue to exist, forever was replaced 
by the notion of a dynamic, expanding universe that seemed to 
have begun a finite time ago. . . " 

Clearly, the expansion of the universe implies that neither spatial scale 
nor time scale remains homogeneous; they both are dynamic. For describ- 
ing the inhomogeneous time, the ideal homogeneous time scale can be used. 
This approach is not new; it has already been used in the theory of relativ- 
ity for describing shortening of time intervals. This means that in fact two 
time scales are considered simultaneously: the ideal, equably flowing homo- 
geneous time, and the cosmic (inhomogeneous) time. The change of scale of 
the cosmic time is described using the homogeneous time scale as a reference 
scale. In other words, the homogeneous time scale is just an ideal notion, 
which is necessary for developing mathematical models describing inhomoge- 
neous cosmic time and its change. In this respect we can, without discussing 
other views on this subject, recall the remark made by A. Daigneault and 
A. Sangalli in their essay [8] about I. E. Segal and his two-time cosmology 
( "chronometric cosmology" , or CC) [25] - note that "perhaps!" : 

"According to CC, Einstein's model is the correct one to under- 
stand the universe as a whole (i.e., global space-time), except 
that there are two kinds of time: a cosmic or Einstein's time t, 
and a local or Minkowski's time xq, which is (perhaps!) the time 
measured by existing techniques. [. . . ] Simply put, Einstein's 
cosmic time is the "real" one, whereas Minkowski's time is only 
an approximation of t." 

So, the ideal model of equably flowing homogeneous time can be considered 
as a rough approximation of the cosmic time. 
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9 Geometric and physical interpretation 
of the Volterra convolution integral 

It should be mentioned that we can also provide a geometric and physical 
interpretation for more general integrals. 

The Riemann-Liouville fractional integral is a particular case of convo- 
lution integrals of the Volterra type: 

K*f(t) = J f(T)k{t-r)dr (24) 


Assuming that k(t) = K'(t), we can write this integral in the form 

/ 

K*f(t) = J f(r)dq t (r), (25) 
o 

q t (r)=K(t)-K(t-r). (26) 

The geometric and physical interpretation of the Volterra convolution in- 
tegral is then similar to the suggested interpretations for fractional integrals. 
The function %(t) determines the changing shape of the "live fence" (in the 
case of the geometric interpretation, see Figs. 1 and 2) and the relationship 
between the individual time and the cosmic time of a moving object (in the 
case of the physical interpretation). 
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